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In this paper, we analyze the feasibility of indoor broadband service provisioning using secondary spectrum 
access to the 960-1215 MHz band, primarily allocated to the distance measuring equipment (DME) system for 
aeronautical navigation. We propose a practical secondary sharing scheme customized to the characteristics of 
the DME. Since the primary system performs a safety-of-life functionality, protection from harmful interference 



becomes extremely critical. The proposed scheme controls aggregate interference by imposing an individual 
interference threshold on the secondary users. We examine the feasibility of large scale secondary access in 
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• terms of the transmission probability of the secondary users that keeps the probability of harmful interference 



below a given limit. Uncertainties in the estimation of propagation loss and DME location affect the feasibility 
of the secondary access. Numerical results show that large number of secondary users are able to operate in 
adjacent DME channels without harming the primary system even with limited accuracy in the estimation of 
the propagation loss. 

Index Terms 

Secondary spectrum access, distance measuring equipment, aggregate interference. 

I. Introduction 

It is generally believed that spectrum shortage is caused by inefficient spectrum utilization under 
the existing regulatory and licensing process that only allows static spectrum allocation. Secondary 
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spectrum access has emerged as a promising solution to relieve the apparent spectrum shortage [Q]. 
In spite of extensive theoretical research on the field of cognitive radio and dynamic spectrum access, 
the practical value of the secondary access has not been fully investigated. Most of efforts to assess 
the real-life benefit of the secondary spectrum have thus far focused on the digital TV broadcasting 
bands, namely TV white spaces 0, 0). Substantial portion of useful spectrum is primarily allocated 
to various systems such as radar and aeronautical navigation, whose secondary access feasibility is 
mostly unexplored 0]. 

This work focuses on the 960-1215 MHz band which is allocated to aeronautical systems. In 
particular, this frequency band is mainly occupied by distance measuring equipment (DME). Secondary 
access to the 960-1215 MHz band was first studied in our previous work [4J. As a first step, we 
investigated the minimum requirements for the secondary users under the ideal assumption that the 
secondary users have accurate knowledge of propagation loss to the DME receivers. We observed that 
the secondary usage would be widely available under this particular assumption. However, it is obvious 
that the requirements to the secondary users will become more stringent if there are uncertainties in 
the propagation information. In practice, it is difficult to have perfect knowledge of the propagation 
to the DME system. Thus, it is needed to study the feasibility of secondary access under practical 
assumptions. 

In this paper, we investigate the feasibility of secondary spectrum sharing with the DME system. 
To our best knowledge, it is the first attempt to examine the practical usefulness of 960-1215 MHz 
with regard to the secondary access. Our contribution can be detailed as follows. First, we propose 
practical methods by which the secondary users discover opportunities and share the spectrum. They are 
customized to the characteristics of the primary user, i.e. DME receivers, and based on geo-location 
database and spectrum sensing. Second, we identify the major sources of uncertainties that cause 
inaccurate estimation of propagation loss to the DME, and analyze the impact of the uncertainties by 
employing mathematical aggregate interference models in fl5), J6]|. 

We consider massive deployment of secondary users that provide high-speed indoor broadband, 
e.g. WiFi and HeNB. Such a large scale secondary access is deemed feasible if the practical sharing 
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methods enable the secondary users to maintain an acceptable transmission probability. Since our 
analysis accompanies the uncertainties in the propagation loss estimation, we focus on the following 
research questions: 

• Is the massive secondary access feasible in 960-1215 MHz band? 

• What is the impact of the uncertainties on the feasibility of secondary access? 

The rest of the paper is organized as follows: the system model, primary and secondary systems 
characteristics are described in Section [I] The proposed secondary access scheme and the mathematical 
models of the aggregate interference for ground transponder and airborne interrogator are introduced 
in Section [TIT] and in Section [IV] respectively. In Section |Vj we present and discuss our numerical 
results. Finally, main conclusions of this work and remaining issues for future studies are given in 
Section ED 

II. System Model 

A. DME as the primary system 

DME is used for measuring the distance between an aircraft and a ground station. The airborne 
equipment (interrogator) sends short Gaussian pulses down to earth, and the ground station (transpon- 
der) responds on a frequency of ±63 MHz from the interrogation frequency. The interrogator can 
calculate the slant distance based on the round trip delay of the signal. The pulses are burst more than 
100 times per second by the interrogator and 2500 times by the transponder. Their transmission power 
reach up to 300 W for the interrogator and up to 2 kW for the transponder. The channel bandwidth 
of DME is 1 MHz, i.e. there are 252 channels in total. More detailed operation of DME can be found 
in BH and references therein. 

We consider that the DME receiver can tolerate a maximum interference power of A t h r , which cor- 
responds to — 119dBm/MHz and — llldBm/MHz for the transponder and interrogator, respectively [|4]|. 
The received interference is considered harmful if it exceeds A thr . Because the DME system performs 
a safety-of-life functionality, protection from harmful interference becomes extremely critical. Due to 
the high sensitivity of DME receivers, we need to control the aggregate interference over a large area, 
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which is the major challenge for the secondary access to this spectrum. The aggregate interference 
(I a ) is regulated as follows: 

Pr[J a > A thr \ < f3 PU (1) 

where (3 PU is the maximum permissible probability of harmful interference at the primary receiver. 
The nature of DME operation requires (3 PU to be extremely small. A reasonable range of (3 PU has 
not been discussed well in the literature. We adopt a value used for air traffic control radar in 2.7-2.9 
GHz, i.e. = 0.001% 0. 

The interference from the DME device to the secondary receiver is, on the contrary, negligible, since 
the DME generates only short pulses. Although the DME pair exchanges the pulses frequently, the 
overall channel utilization remains below 1%. Secondary receivers might be saturated if they receive 
excessively strong DME pulses. Let I sat be the saturation point of the secondary receiver. Then, the 
following condition should be satisfied: 

Pt[I PV > I sat ] < PSU (2) 

where fisu is the maximum saturation probability and I P u is the received primary pulse power. We 
adopt a value of psu — 2% and I sat = — 30dBm which is a typical saturation level of low noise 
amplifier (LNA) in WiFi receivers [[8]. With the adopted values for A thr , a simple link budget analysis 
indicates that (OQ) is the limiting constraint even before taking the effect of multiple secondary users 
into account. Therefore, we will focus on the protection of the primary user in the remainder of the 
paper. 

B. Indoor Broadband as secondary system 

Let us consider massive scale deployment of indoor access points and mobiles for high capacity 
broadband services over a large area. They utilize the spectrum allocated to the DME under the principle 
of spectrum interweave [0. The secondary users are assumed to be spatially distributed according to a 
homogeneous Poisson point process in a two dimensional plane !R 2 . The primary receiver is located at 
the center of the circular region limited by two radii r a and R, which are the minimum and maximum 
distances from the primary receiver, respectively. 



5 



Each secondary user decides whether it can access a particular DME channel or not by estimating 
the interference it will generate to the primary user. Let I t hr denote the interference threshold imposed 
on the individual secondary users. The value of I thr is given to the secondary users by a central 
spectrum manager. This ensures that each secondary users makes its own decision without interacting 
with the others. The interference from a secondary user i is given by 



where & is the interference that the primary user would receive if an arbitrary secondary user were to 
transmit, and & is the estimate of & by the secondary user i. Note that & = £j only when the secondary 
user has the perfect knowledge of the propagation loss. Considering that there are N secondary users 
around the primary user, the aggregate interference is 



A. Secondary access scheme 

The ground transponder is placed at a fixed location and frequently bursts short pulses to the airborne 
interrogators. Thus, it is possible for the secondary user to detect the existence of the transponder via 
spectrum sensing. The additional use of geo-location database enables the secondary users to have 
prior knowledge about the transponder such as the location, operating frequency, and transmission 
power. This will significantly improve the performance of the spectrum sensing since the secondary 
users can have a good expectation about to signal to detect. Given the high transmission power of the 
transponder, we assume that the spectrum sensing is reliable enough to ignore missed detection and 
false alarm. 

Fig. Q] depicts the proposed opportunity detection mechanism. Notice that the secondary users detect 
the transponder on the reply (sensing) frequency, while the interference is given on the interrogation 
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where N t is the set of transmitting secondary users. 



III. Sharing with the ground transponder 



(interfering) frequency. In both channels, propagation losses between the DME transponder and the 
secondary user consist of the distance-based path loss (L) and fading^ (X and Y). Although it is 
reasonable to assume that the secondary users accurately estimate the propagation loss of sensing 
channel (S = L+X), it does not necessarily mean that the estimation of interfering channel (T = L+Y) 
is also accurate. With the frequency offset of 63 MHz between the sensing and interfering channels, 
the shadowing components are typically highly correlated (p s hadowin g ~ 1) ifTOl . while the multi- 
path fading is uncorrected (pf as t = 0). Therefore, the correlation between the composite fading 
components, p, lies between [0, 1]. The exact value of p depends on the characteristics of different 
propagation environments. Partial correlation between channels does not allow the secondary user to 
perfectly estimate its interference to the primary victim. Then, an uncertainty in the estimation of 
fading component of the propagation loss between the secondary user and the ground transponder still 



remains. 



B. Aggregate interference modeling 

In this section, we model the aggregate interference when there is uncertainty in the fading estima- 
tion. Different levels of uncertainty in fading estimation are represented by a correlation coefficient p. 
We adopt the mathematical frameworks proposed in []4)-[]6] with a slight modification to account for 
the proposed spectrum sharing mechanism. 

Let us consider an arbitrary secondary user i which is distributed according to a homogeneous 
Poisson point process in a circular area of radius R. The path loss between the primary receiver and 
the secondary user i is modeled as g(rj) = Cr^ a where C is a constant and a is the path loss exponent. 
Then, the user % would cause interference ^ to the primary receiver if it were to transmit, which can 
be expressed as 

& = P t eff g(r t )Y t (5) 

where P^^ refers to the effective transmission power of the secondary user including antenna gains and 
bandwidth mismatch. Y { is a random variable modeling the fading effect. It is generally considered that 

'Note that the fading here refers to the combined effect of shadowing and multi-path fading 



the fading consists of shadow fading following a normal distribution in dB scale and multi-path fading 
by which the instantaneous power is varied with an exponential distribution. We use the assumption 
that the composite fading Y { follows a log-normal distribution. It is known that this assumption works 
well when the standard deviation of shadowing is higher than 6dB, i.e. when the shadowing is a 
dominant factor of the composite fading [fTTfl. 

The user i will decide to transmit if £j < I t hr- Note that £j is affected by the fading on the sensing 
channel. That is, 

6 = Pt ff g{ri)Xi (6) 

where X { is modeled as a log-normally distributed random variable whose parameters are same as Y^. 
Therefore, the joint distribution of X,- L and Yi is given by the following bivariate log-normal distribution: 

1 (ln:c) 2 -2p(ln 3: )(ln i) ) + (lni ; ) 2 

fxtx&y) = 2 n f e ^ IT ^ T) (7) 

Zitxya^yl — p z 

where p is the correlation coefficient of Xi and Yf. 

^Var (In X~) Var (In Y~) 

We consider that the composite fading components, Xi and Y { , will be partially correlated (0 < 
p < 1). The exact value of p depends on propagation environments. Note that full correlation (p = 1) 
represents an ideal case that the secondary user has an accurate knowledge of interference. On the 
opposite, zero correlation (p = 0) stands for a pessimistic assumption that the fading is completely 
unknown to the secondary user. For simplicity and mathematical tractability, we have adopted the 
assumption that secondary users in the whole area of study are affected by a homogeneous fading 
distribution. The feasibility of secondary access under different assumptions, ranging from ideal to 
pessimistic, will be shown and discussed in Section |V] 

The aggregate interference I a can be expressed as: 

I a = P t eff Cj2 r i a Y- (9) 




Hereafter, we omit the index of secondary user i, which is chosen in an arbitrary manner, unless 
necessary. By applying the Campbell's theorem, the characteristic function of Jjv t is as follows: 

^iN t ti w ) = ex P ( ~ 2nX [ [ I [1 - exp(jwyr- a )} 

\ JX JY J To 

x \ 0tI ; hr] (r~ a x)f x ,Y(%, y)rdrdydx^ . 

where j = \/—\ and I t h r = hhr/{Pt C). The activity of the secondary users is represented by 
lmj^ Ar~ a x), which is a Bernoulli random variable. The indicator function is defined as: 

{1, if a < z < b 
(11) 
0, otherwise 

where the value one of the Bernoulli variable denotes that the secondary user is able to transmit. We 
use (fTOl) to derive exact expressions for the n th cumulant of the aggregate interference in a limited 
circular region [r , R]. We consider the case where there is a partial correlation between the two fading 

effects affecting the sensing and interfering channels, X and Y . 
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where, 



_ ln(rg/ tfer ) - piny 
~ ^ ' 
_ \n(R a It hr ) - piny 

L, — = . 

For the special cases of full correlation (p = 1) and zero correlation (p = 0), the closed-form 
expressions of cumulants can be found in [5 J and [0, respectively. Using the cumulant of shown 
in (PT21) . we can obtain the n th cumulant of the aggregate interference I a as follows: 

k Ia (n) = (Pt eff C) n k lNt (n). (13) 



The probability density function (pdf) of I a can be approximated with a known distribution by 
moment-matching method. In Q, J6l, shifted log-normal and truncated-stable distributions are em- 
ployed to address the skewness of the aggregate interference. In our model, the strong interferers are 
effectively removed due to the stringent threshold in (0. Therefore, simple log-normal distribution 
sufficiently describes I a . The pdf of I a can be approximated with the first and second order cumulants 
of I a obtained by (fl3i 



fia(y) = — 7^= exp 



(14) 



where 



k Ia (l) = ex P [ f i Ia + al/2], (15) 
k Ia (2) = exp(4 - 1) exp(2/i /a + a 2 J. (16) 

IV. Sharing with the airborne interrogator 
A. Secondary access scheme 

Airborne interrogators are equipped in the airplanes, which are moving with a high speed. Therefore, 
it is not reasonable to assume a reliable detection of the interrogator via spectrum sensing. Instead, 
we assume that the secondary users are connected to a real-time database where the locations of the 
airplanes are provided. A living example of such a real-time aircraft location map can be found in 
lfl2l . Currently, the database information is updated every 20-60 seconds and has a limited coverage, 
which means that some airplanes (mostly small ones) do not appear in the map. However, we expect 
that an official database in a future will be able to provide a reliable information since it will be 
maintained by national authorities. 

Due to the update delay in the database, the secondary user could potentially experience uncertainty 
or imperfect information on the location of the airborne interrogator which is changing rapidly. Based 
on the update delay and the speed of the airplane, we introduce the notion of error region, inside 
which secondary users will assume the worst case scenario that the sky is full of airplanes as shown 
in Fig. El Outside the error region, secondary users will assume that the primary receiver is located 
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at the closest border of the error region. Let t u be the time of update delay and v be the speed of 
the airplane. Then, the radius of the error region is given by t u v. For example, The t u of one minute 
corresponds to the error region of 15 km radius assuming v = 900 km/h. 

B. Aggregate interference modeling 

For the case of the airborne interrogator, free-space propagation model between the secondary users 
and the primary receiver is assumed. This means that fading effect is not taken into account. We adopt 
this assumption in order to account for the worst case scenario where there exists line-of-sight path 
between every secondary user and the primary user. 

Similar to the ground transponder case, we assume that N secondary users are distributed according 
to a homogeneous Poisson point process in a circular area of radius R. The primary victim is assumed 
to be located at the center with a height of h from the ground. Since the fading effect is not considered, 
applying individual threshold I thr will result in a circular exclusion region where secondary users are 
not allowed to transmit. The radius of exclusion region is denoted by r Q . 

Let r thr be the exclusion radius under the assumption that the secondary users know the exact 
location of the primary victim. In the presence of the update delay, each secondary user has to make 
a conservative decision that the airplane is at the closest possible location. It effectively increases the 
exclusion radius by t u v. However, if the exclusion region is not needed in the first place (r thr = 0), the 
uncertainty in the primary user location does not make any impact on the feasibility of the secondary 
users. Thus, r Q is given by 



Let li be the distance from an arbitrary secondary user % to the primary receiver. Then, U = y/h 2 + 
and the path loss g(li) is given by Cl7 L a '. Then, the aggregate interference I a is 

N t 

I a = P t eff Cj2^ a (18) 

i=i 





otherwise. 



(17) 
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where N t is the number of secondary users that are allowed to transmit, i.e. located at the outside 
of the exclusion region. Similar to Section IIII-BL we apply the Campbell's theorem to obtain the 
characteristic function of Jjy t . Then, we derive exact expressions for the n th cumulant of in a 
limited circular region [r , R\. 

klN ( n ) = J^( B (2-n«)/2 _ A (2-na)/ 2) _ (19) 

' net — 2 

where A = h 2 + R 2 and B = h 2 + r 2 . Since we consider the free-space propagation model (a = 2), 
we employ l'Hopital's rule to calculate the first order cumulant (kj N (1)). Using the cumulant of I Nt , 
we can obtain the n th cumulant of the aggregate interference I a as it is shown in (fl"3l) . Note that Jj is 
only affected by the distance-based path loss. Thus, I a is well described by the central limit theorem. 
This means that I a can be approximated as a Gaussian distribution with the first two cumulants as the 
mean and variance. 

V. Numerical Results 

The parameters used for our numerical experiments are described in Table HI For the case of ground 
transponder, we model the propagation loss between the primary victim and the secondary user using 
Hata model for suburban area. Instead, for airborne interrogator we employ free-space propagation 
loss. For the transponder, we investigate the impact of p on the requirement and feasibility of secondary 
access in terms of the individual interference threshold I t h r an d the transmission probability of the 
secondary user i, Pr(£j < I t hr), at a given r^. For the interrogator, we analyze the effect of the update 
delay on the requirements of secondary users. The feasibility of secondary access is given in terms of 
the exclusion region size r Q imposed on the secondary users. 

In both cases, we provide results for co-channel usage and as well as adjacent channel usage. 
We apply DME selectivity mask given in 1(131 to determine the adjacent channel rejection (ACR) 
characteristics. This means that the condition (H) is changed to Pr[/ a > (A t h r + ACR)] < (3pu 
when we evaluate the adjacent channel usage. The values of ACR will vary according the frequency 
separation. For instance, it is between 60dB and 70dB for channels with a frequency separation of 2 
MHz. We assume that this applies as well to the channels of more frequency separation. To account 
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for interference aggregation in the spectral domain, we apply a fixed margin of 3dB and lOdB for 
co-channel and adjacent channel, respectively. 

For the case of ground transponder, the cumulative distribution function (CDF) of I a calculated from 
(fl4)) for different values of p is shown in Fig. [3] A good agreement between analytical CDF of I a 
and the simulation results is verified when p > 0. When the fading is unknown to the secondary user 
(p = 0), analytical CDF matches the tails of the simulation-based CDF of I a . Since we are working 
with (3pu = 0.001%, it is still possible to employ the log-normal approximation of the probability 
distribution of I a to analyze the impact of fading uncertainty on the feasibility of secondary access. 

The individual interference threshold I t hr required to satisfy (H) is given in Fig. |4] It shows the 
impact of p on the required I t hr for accessing a co-channel and an adjacent channel with ACR of 
60dB. We observe that the margins to cope with the uncertainty for different values of p does not change 
much when the density of secondary users per km 2 (Xsu) increases. However, the uncertainty margin 
significantly varies for co-channel and adjacent use cases, i.e. for different A t h. r values. Considering 
that secondary users transmit in an adjacent channel with ACR of 60dB, it is observed in Fig. [5] that 
the impact of fading uncertainty is critical for high-power secondary users (above lOdBm). However, 
the operation of a dense secondary network for indoor coverage is feasible even if the secondary users 
cannot accurately estimate the propagation loss. 

Now, let us consider the airborne interrogator as the primary victim. Recall that the update delay 
of 5 minutes can lead to the error region of 75 km radius, which is almost equivalent to not having 
the database. The exclusion region needed to satisfy (ff) is shown in Fig. [6] The impact of the update 
delay is significant only when ACR is lower than 50dB. When ACR is 60dB, no exclusion region is 
required even if long update delay is experienced in the communication between the secondary user 
and the real-time database. Fig. |7] shows the combination of secondary users density and transmission 
power that do not require fast database update, i.e. no requirement for exclusion region. The figure 
indicates that dense secondary network accessing adjacent channels is feasible when the transmission 
power is about OdBm even if no information on the location of the primary victim is provided. 
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VI. Conclusions and Future Work 

We analyzed the feasibility of large scale indoor broadband secondary access to the 960-1215 
MHz spectrum when uncertainties on the fading and the location of the primary receiver are present. 
Cumulant-based approximations have been employed to derive the probability distribution of the 
aggregate interference in the presence of uncertainties. The main contributions in this paper are twofold: 

• We proposed a practical secondary sharing scheme considering the characteristics of different 
primary receivers (DME ground transponder and airborne interrogator). Then, we identified un- 
certainties in the estimation of propagation loss incurring from the proposed sharing scheme. 

• The feasibility of large scale secondary access has been evaluated in terms of the number of 
secondary users which are able to operate with an acceptable transmission probability and the 
exclusion region size imposed on the secondary users. 

We conclude that massive indoor secondary access to adjacent channels (ACR > 60dB) is feasible 
even if secondary users are not capable of accurately estimating the propagation loss nor have accurate 
knowledge of the location of the airborne interrogator. Numerical results show that dense secondary 
users (\su > 1000/km 2 for ground transponder and \ S u > 100/km 2 for airborne interrogator) can 
have access to adjacent channels with a high transmission probability (> 90%) or small exclusion 
region size. 

Since the indoor secondary use of 960-1215 MHz spectrum is identified feasible, the capacity 
analysis of the secondary system taking self-interference and power control into account remains as 
an interesting future work. Location-dependent availability of the secondary access and its economic 
value are also to be investigated. 
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TABLE I 

Parameters used for numerical experiments 



Parameters for ground transponder 



primary user transmission power 


60dBm/MHz 


path loss constant (C) 


4.5 x lCT 13 


path loss exponent (a) 


3.5 


Fading standard deviation (<r^) 


lOdB 


height of the transponder 


30 m 


Parameters for airborne interrogator 


primary user transmission power 


55dBm/MHz 


path loss constant (C) 


5.7 x 1CT 10 


path loss exponent (a) 


2.0 


height of the interrogator (h) 


1 km 


Common parameters 


radius of interference aggregation (R) 


200 km 


building penetration loss 


lOdB 


DME antenna gain 


5.4dBi 


secondary user antenna gain 


OdBi 


secondary user transmission power 


ldBm/MHz 


secondary user height 


1.5 m 
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Fig. 1. Uncertainty in secondary sharing scenario with ground transponder as primary victim 




Fig. 2. Uncertainty in secondary sharing scenario with airborne interrogator as primary victim 
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Fig. 4. I t hr as a function of secondary network density (Xsu) for different correlation coefficients(p); the primary receiver is the DME 
ground transponder 
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Adjacent Channel, ACR = 60dB 




^/Km 2 



Fig. 5. Impact of different correlation coefficients (p) on the feasibility of the secondary access when Pr(£j < Ithr) > 90% at 
rj — 5km, the primary receiver is the DME ground transponder 




Fig. 6. Exclusion radius as a function of adjacent channel rejection of DME for different update delays when \su = 20/km 2 , the 
primary receiver is the DME airborne interrogator 



19 



Adjacent Channel, ACR = 60dB 

30 | 1 1 1 , 
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Fig. 7. Maximum secondary user transmission power for a given Xsu when no exclusion region is needed (r = 0km), the primary 
receiver is the DME airborne interrogator 



